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Summary 



The primordial abundances of deuterium and of helium-3, produced during big 
bang nucleosynthesis, depend sensitively on the baryon density. Thus, the ob- 
served abundances of D and 3 He may provide useful "baryometers" provided the 
evolution from primordial to present (or, presolar nebula) abundances is under- 
stood. Inevitably, the derivation of primordial from observed abundances requires 
the intervention of a model for galactic evolution and, so, the inferred primordial 
abundances are, necessarily, model dependent. Here, an analytic framework for 
the evolution of D and 3 He is presented which is "generic" in the sense that it 
should describe the results of any specific galactic evolution model. The "effective 
3 He survival fraction", T^, is the one free parameter which is model specific. So- 
lar system and interstellar data are used to infer upper and lower bounds to the 
primordial deuterium mass fraction (X^p) as a function of T3 and, these bounds 
are used to constrain the present baryon-to-photon ratio (rj) and baryon density 
For r 3 > 1/4 it is found that (from D and 3 He alone): 3.1 < rji < 9.0; 
0.045 < n B h 2 50 < 0.133 (where H Q = 50h 50 kms" 1 Mpc" 1 ). 
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Introduction 

The light nuclides D, 3 He, 4 He and 7 L i are produced in astrophysically in- 
teresting abundances during big bang nucleosynthesis (BBN) and, hence, their 
primordial abundances provide quantitative tests of the standard hot big bang 
model. However, primordial abundances are not observed but, rather, derived 
from observational data and - to a greater or lesser extent - models for galactic 
chemical evolution. Since 4 He and 7 Li are observed in metal-poor objects (extra- 
galactic HII regions for the former, PopII stars for the latter), the evolutionary 
corrections are minimized. Unfortunately, this is not the case for D and 3 He where 
observations are, for the most part, limited to the solar system and the present 
interstellar medium (ISM) . The inevitable intervention of galactic evolution mod- 
els results in model dependent uncertainties for the primordial abundances of D 
and 3 He derived from even the most accurate observational data. Indeed, for this 
reason the most fruitful approach has been to use the observational data to bound 
the primordial abundances. For example, since deuterium is only destroyed dur- 
ing the course of galactic evolution, the primordial abundance must exceed the 
observed abundance (X2P > -V20BS where X20BS = -^2© or X2ism)- It has, how- 
ever, proved harder to use the observational data to bound X2P from above (Yang 
et al. 1984, YTSSO). YTSSO and Dearborn, Schramm & Steigman (1986, DSS) 
first noted that since D burns to 3 He and some 3 He survives stellar evolution, the 
observed abundances of D and 3 He may be used to bound, from above, the primor- 
dial abundance of D + 3 He . The upper bound to the primordial abundance of D 
+ 3 He derived by YTSSO, and utilized in many subsequent analyses, depends on 
the 3 He survival fraction (73 which DSS computed from models for the evolution 
of stars of different masses. However, unlike the model-independent lower bound 
(X 2 p > ^20Bs), this upper bound to D + 3 He is still model dependent in the 
sense that the "effective" value of (73, hereafter indicated by 1?3 to distinguish it 
from the DSS g3, depends not only on stellar structure/evolution but, also, on 
the past star formation history, the IMF (and its evolution), infall and outflow, 
etc. In short, since the "true" value of T 3 is model dependent, so too is the upper 
bound to the primordial abundance of deuterium derived from observations of D 
and 3 He. Indeed, using very much the same stellar data, Steigman & Tosi (1992, 
ST92) find deuterium destruction factors < 2-3 while Vangioni-Flam & Audouze 
(1988) and Vangioni-Flam, Olive & Prantzos (1994) find destruction factors as 
large as ~ 5 — 10 are allowed and, the more extreme models of Vangioni-Flam & 
Casse (1994) permit even more D destruction. We note that these latter mod- 
els have been designed to account for the D data and they may not necessarily 
be consistent with other observational constraints on the abundances, abundance 
ratios, and gradients of the other elements. Indeed, Edmunds (1994) argues on the 
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basis of "Simple Model" predictions, that they may not be consistent with the 
age-metallicity relation and the G-dwarf metallicity distribution. 

Despite the above-cited model dependent uncertainties, it is important to 
pursue the primordial abundance of deuterium. The main reason is that, in 
the context of BBN, X2P (the primordial deuterium mass fraction) is a sensi- 
tive "baryometer" . At higher nucleon (baryon) densities D burns faster (to 3 H, 
3 He, 4 He) and less deuterium emerges from BBN. X2P is a steeply decreasing 
function of the nucleon-to-photon ratio ?y[?7io = 10 10 (iV/7)]; for 2 < rj\Q < 6, 
X2P ~ Viq 7 - Thus, an upper (lower) bound to X2P will lead to a lower (upper) 
bound to r\ and, hence, to bounds on the current universal density of baryons 
(Qb = PbI Vcrit = 0.015?7io/i5 O 2 where H = 50/i 50 km s _1 Mpc -1 ). Furthermore, 
a lower bound to r\ (derived from an upper bound to X 2 p), when combined with an 
upper bound to the primordial 4 He mass fraction, can be used to bound the "effec- 
tive" number of equivalent light neutrino degrees of freedom (Steigman, Schramm 
& Gunn 1977) and, hence, to constrain "new" physics (beyond the standard model 
of the strong and electroweak interactions). Thus, rather than ignore primordial 
deuterium as a potential baryometer, it is surely more worthwhile to try to under- 
stand the nature and extent of the model dependent uncertainties in its inferred 
primordial abundance. To this end, it is the goal here to derive simple analytical 
relations between the primordial and observed abundances of D and 3 He which 
are general in the sense that they do not rely on the special details of any specific 
chemical evolution model. The equations to be derived in the next section are 
"generic" but, they depend on one model dependent parameter, 1?3, the "effective 
survival fraction of 3 He". Different galactic evolution models will yield different 
values of T 3 (as a function of time and of location in the Galaxy). Highlighting the 
physical reasons for the various values of T 3 may aid in identifying a "reasonable" 
range for Ts and, hence, lead to nearly model independent bounds on X2P, 1710 
and Vtp- 

Evolution of Deuterium and Helium-3 

At any time (t) consider a "representative" volume (V) of the ISM which, 
initially, contained an amount of gas Mj. In general, the gas mass in V will change 
with time due to infall, outflow and the sequestering (even if temporary) of gas in 
stars. At any time during the evolution of the Galaxy there will be a fraction, f(t), 
of the gas in the ISM (in V) which has never been through stars. The remaining 
fraction, 1 — f(t), has been cycled through stars - one or more times - and returned 
to the gas. Thus, at time t in volume V, f(t)M(t) has never been through stars 
and (1 — f(t))M(t) has been through one or more generations of stars. Further, 
of the amount of gas in V at t, M(t), there were initially (prior to any evolution) 
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N 2i deuterium nuclei where, 

N 2i = (X 2i /2)Mi/M N . (1) 

In (1), X 2 i is the initial (primordial) D mass fraction and Mn is the nucleon 
mass. Since any D cycled through stars is destroyed, at time t in volume V only 
N 2 (t) D-nuclei remain. 

N 2 (t) = f(t)(X 2i /2)(M(t))/M N . (2) 

Since X 2 {t) = 2N 2 (t)MN /M(£), the D-survival factor (f 2 ) and the virgin gas 
fraction (/) are identical. 

f 2 (t)=X 2 (t)/X 2i = f(t). (3) 
Since deuterium is only destroyed, f2(t) < 1. 

Next, consider the evolution of 3 He. 3 He is destroyed in the hot interiors 
of all stars but preserved in the cooler outer layers (Iben 1967; Rood 1972). For 
cooler low mass stars (< 2.5M ) hydrogen burning results in the net production 
of 3 He (Iben 1967; Rood 1972; Rood, Steigman & Tinsley 1976; DSS). Although 
net stellar production of 3 He is potentially large, especially in the epoch between 
the birth of the solar system and the present, when the most effective 3 He pro- 
ducers start to contribute to the ISM enrichment, it is also very uncertain (Rood, 
Steigman & Tinsley 1976; DSS; Galli et al. 1994, Hogan 1994, Tosi et al. 1994). 
Thus, in the spirit of deriving "conservative" upper bounds to the primordial 
abundances of D and 3 He, newly synthesized 3 He (as distinct from prestellar D 
burned to 3 He) is ignored in the following. This will lead to a lower bound to 
X 3 (t), the 3 He mass fraction at any time t. The quantitative effect of this neglect 
is reduced if presolar nebular abundances are used rather than current interstellar 
values. Now, since a fraction f(t) of the 3 He nuclei which were originally in M(t) 
have never been in stars and a fraction 1 — f(t) have been cycled through stars, 
for X 3 (t) = 3N 3 (t)M N /M(t) it follows that, 

X 3 (t) > X 3l f(t) + (1 - f(t))T 3 (t) [X 3l + 3X 2l /2] . (4) 

In (4), T 3 (t) is the "effective survival fraction" of 3 He. 1?3 is model and time 
dependent since it depends not only on the stellar models but, on the IMF and 
on the history of star formation (e.g., which stars have returned their material to 
the ISM by time t). In addition, although in material going through stars for the 
first time the prestellar 3 He abundance is the sum of the D and 3 He primordial 
abundances (D is immediately burned to 3 He), any material going through stars 
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for a second (or more) time will have had its initial deuterium already destroyed; 
this, too, modifies r 3 from the simple DSS values g 3 and, thus, r 3 depends on 
how much of the gas that has been through stars at least once has been through 
subsequent generations. 

The 3 He evolution factor /3(f) is, 

f 3 (t) = X 3 {t)/X 3l > f(t) + (1 - f(t))T s (t)(y 23i /y Si ). (5) 

In (5), y 23 i/y 3 i is the initial ratio (by number) of D+ 3 He to 3 He. By combining 
equations (3) and (5) we may eliminate fit) and directly relate X 2 i and X 3 i to 
X 2 (t) and X 3 (t). However, before doing so, it is instructive to explore those 
conditions under which X 3 will increase from its initial value (/3 > 1). Eq. (5) 
may be rewritten as, 

fs(t) > 1 + (1 - /(f)) [TaWd/Mi/vsO - 1] • (6) 

Since 1 — f(t) > 0, the 3 He mass fraction will increase if T 3 (t) > y 3 i/y 23 i- As 
an illustration, if T 3 > 1/4 (DSS), (y 23 /y 3 )p > 4 is required; for standard BBN 
(Walker et al. 1991, WSSOK) this is satisfied for r] W < 4.5. Thus, if r 30 > 1/4 
and 7710 < 4.5, it is predicted that X 3& > X 3 p(> 2.9 x 10~ 5 ); this is (as will 
be seen) consistent with solar system data (Geiss 1993). Alternatively the ST92 
models 1 suggest that r 30 > 1/2 which implies X 3& > X 3 p for all 7710 < 10. 

The model dependent factor f 2 (t) = f(t) may be eliminated from (5) by the 
use of (3), relating the abundances (by number with respect to hydrogen) of D 
and 3 He initially to those at any later time t. 



V2i < 



X 

x~ 



2/23 + 



1 



-1)2/3 



1 



V23i 



2/23 



(7) 



In (7), X = X(t) is the hydrogen mass fraction and Xi is its initial value. The 
inequality is due to the neglect of stellar synthesis of 3 He and is generic to all 
models of galactic evolution. The only model dependence in (7) is in T 3 (t). It will 
be noticed that eq. (7) is a "mixed" expression in the sense that the last term on 
the right hand side depends on both initial and final abundances. One approach 
(YTSSO) is to neglect this term, strengthening the inequality; this is equivalent 
to assuming zero initial 3 He. For y 3i > 0, 



2/2; < y 2i < 



x 



2/23 + 



1 



- 1 



2/3 



(8) 



1 Due to unfortunate confusion in the DSS definition of g 3 (which is by mass 
and not by number), the survival fractions f 3 published by ST92 were underesti- 
mated on average by a factor 1.2. This implies r 30 R3 T 3I sm ~ 0.60 
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Except for the prefactor (X/Xi), eq. (8) (for y^) is precisely the expression 
derived in YTSSO and used in many subsequent analyses. However, this prefactor 
is required to account properly for the changing hydrogen mass fraction. Since 
X(t) < Xi the "usual" bound to y^i is somewhat overestimated; for example, for 
X P = 0.76 ± 0.02 and X(t) = 0.70 ± 0.02, X{t)/X P = 0.92 ± 0.04. 

The symbiotic relationship between D and 3 He is clear in equations (7) and 
(8) where, for fixed (observed) abundances yi and 2/3, the initial abundances of D 
and 3 He are anticorrelated (e.g., higher y^i forces lower 2/2i)- Although eq. (8) 
provides a valuable BBN-independent bound to pre-galactic deuterium (Steigman 
1994), the last term on the right hand side of (7) should not, in general, be 
neglected. This poses no problem in practice since, for a specific BBN model the 
left and right hand sides of (7) may be compared to find if the inequality is, or 
is not, satisfied. In this manner the observational data can be used, along with 
evolution model-dependent choices of 1?3, to constrain the models/parameters of 



In subsequent comparisons, the solar system data of Geiss (1993) will be 
adopted. The observed quantities, derived from solar wind, lunar and meteoritic 
data are the ratios (by number) to 4 He of presolar D ± 3 He and 3 He (1/23/2/4 and 
2/3/2/4)- With such data it is convenient to rewrite (7) and (8) in the following 
forms, 



In (9) and (10), Y Q is the solar 4 He mass fraction. 

Solar System Abundances 

In a very valuable recent review, Geiss (1993) has reanalysed the solar sys- 
tem data crucial to determining the presolar nebular abundances of D and 3 He. 
Although Geiss' central values are virtually identical to those used in the past 
(e.g., Boesgaard & Steigman 1985), his error bars are more generous. Here, with 
a little care towards error propagation, the Geiss (1994) analysis is adopted. 



BBN. 




(2/23/2/4) 



= 4.09 ±0.83 x 10 



-4 



(11a) 







(2/3/2/4)© = 1-52 ±0.30x10 



-4 



(116) 
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The adopted value of the solar 4 He abundance, with generous allowance for un- 
certainty, is 

Y Q = 0.28 ±0.02 (y 4Q = 0.10 ±0.01). (12) 
For Z Q 0.02, X Q = 0.70 ± 0.02 and the above values may be used to derive: 



(Y Q /2)(y 23 /y 4 ) Q 


= 5.73 ± 1.23 x 10" 


-5 


(13a) 


(Y Q /2)(y 3 /y 4 ) Q 


= 2.13 ±0.45 x 10" 


-5 


(136) 




= 3.60 ± 1.26 x 10" 


-5 


(13c) 


X30 


= 3.19 ±0.67 x 10" 


-5 


(13d) 


2/23© 


= 4.09 ±0.92 x 10" 


-5 


(13e) 


2/3© 


= 1.52 ±0.34 x 10" 


-5 


(13/) 


2/2© 


= 2.57 ±0.92 x 10" 


-5 


(1327) 



Bounds To Primordial Deuterium 

Armed with solar system abundances, equations (3) and (9) (or (10)) may be 
enlisted to provide upper and lower bounds to the pregalactic (primordial) abun- 
dance of deuterium. If equation (9) is used, y 3 i/y 23 i will be taken from standard 
BBN (e.g., WSSOK). As already emphasized, the only galactic evolution model 
dependence is contained in r 30 , the effective 3 He survival fraction in the solar 
neighbourhood at the time of the formation of the solar system. For quantitative 
comparisons Ts & = 1/4 and T 3 q = 1/2 will be adopted; in general, the larger r3©, 
the smaller the upper bound to X 2 p. 

Since, at the time of formation of the solar system, at least some of the 
material then present in the local ISM had been cycled through stars, /q < 1. 
Thus, X 2 p > X 2Q ; this differs from the often quoted constraint y 2 p > y 2 Q by the 
factor Xq/Xp ^ 0.9. At the 2-sigma (~ 95% confidence) level, it follows from 
(13c) that, 

X 2P > X 2Q > 1.08 x 10" 5 . (14) 

Because of the relatively large uncertainty in X 2Q , this only provides a very weak 
constraint on X 2 p. For a more restrictive lower bound to primordial deuterium, the 
ISM D/H ratio (Linsky et al. 1993) is more useful. The interstellar observations 
from Copernicus and IUE (see, e.g., Boesgaard & Steigman 1985) and HST (Linsky 
et al. 1993) extend over two orders of magnitude in H column density and are 
consistent with 

z/ 2ISM = 1.6±0.2xlO- 5 . (15) 
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If X lSM » X, 



© 



0.70 ± 0.02, then at 2-sigma, 



X 2P > X 2 isM > 1.7 x 10 



-5 



(16) 



For standard BBN (WSSOK) this provides an upper bound to the nucleon density, 



where 77 is the ratio of nucleons to photons at present and 7710 = 10 77. 

As an aside, notice that if the solar system (13g) and interstellar (15) D/H 
ratios are compared, 



This provides a weak (~ 1-sigma) hint that, in the solar neighbourhood, some D 
has been destroyed in the last 4.6 Gyr. 

In considering the upper bound to the primordial deuterium mass fraction 
from equations (9) or (10), it is useful to compare the BBN predicted value A^bbn 
(as a function of 77) with the analytic bounds X 2 uax (eq. 9) and J^max ( ec l- 
10). For a "2a" upper bound, equations (9) and (10) are evaluated with the 2a 
upper bounds to the solar system abundances (see 13a & b: (Yq/2) (2/23/2/4)0 < 
8.2 x 10" 5 , (Yq/2) (t/sM)© < 3.0 x 10" 5 ). The results for X 2MA x and X 2 ° MAX 
are shown in Figure 1 for T 3 = 1/4 & 1/2 respectively. Also shown in Figure 1 
is X 2 bbn] the requirement that X 2 bbn < ^2MAX (or J^max) provides a lower 
bound to r/. 

As anticipated, the neglect of any primordial 3 He (^max) provides a weaker 
constraint than if the BBN predicted relative 3 He/D abundance is allowed for 
(A^2max); for all 2/3P > 0, -X^max > ^2MAX • As 77 increases, the BBN yields of 
D and 3 He decrease but, ( 3 He/D)BBN increases. It is the increasing importance of 
3 He at high 77 which is responsible for the decrease in X 2 max ■ 

As more 3 He survives stellar processing increases), consistency with the 
solar system abundances restricts primordial D to lower values (X 2 max and ^max 
decrease). Thus, for fixed observed abundances, the lower r 30 , the weaker the 
constraint on X 2P . For r 30 > 1/4, X 2P < X$ MAX < 16 x 10" 5 (@ 2a) which, 
for BBN, corresponds to 77^ > 2.5. However, for 7710 > 2.5, 2/3p/2/23P > 0.14 and 
the neglect of primordial 3 He is non-negligible. From Figure 1 it may be seen that 
for r 30 > 1/4, A2BBN < I2MAX (@2a) for X 2P < 11 x 10~ 5 ; this corresponds 
to 77iq IN > 3.1. For the same solar system data, if r 30 increases to > 1/2, ?7{^ IN 
increases to > 4.0. 



7710 < 9.0, 



(17) 



^2©/^2ISM ~ 2/2©/2/2ISM = 1-6 ± 0.6. 



(18) 
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In Figure 2, X 2 bbn is compared to the 2a upper bound to X 2 max for T 3Q = 
1/4 and, to the 2a lower bound to X 2 mw from the interstellar D/H ratio (eq. 
16). From D and 3 He alone, r\ is then restricted to the range 3.1 < 7710 < 9.0. 
Observational data on 4 He and 7 Li (e.g. WSSOK; Olive & Steigman 1994) further 
constrain the upper bound to 77(7710 < 3.9). 

Earlier, it had been noted that, provided T 3 is not too small nor y 3 p/y 23 p too 
large, it is likely that the 3 He mass fraction will increase from its primordial value 
(even neglecting net stellar production of 3 He). Thus, for F 3 @ > 1/4 and 7710 < 4.5, 
it was expected that X 3 q > X3bbn- This upper bound to primordial 3 He will also 
lead to a lower bound to r\ which, although weaker than that above, is still of some 
interest. From (13d), @2cr, X 3P < X 3& < 4.5 x 10~ 3 , which corresponds to 77 10 > 
2.0. Since this bound is independent of the previous deuterium constraint it may 
be used to bound primordial/pregalactic D: for 7710 > 2.0, (D/H)p < 1.6 x 10 -4 
(and, for the standard model of 3 light neutrino species, Yp > 0.236). This bound 
is of some interest considering the claim of a possible detection of deuterium in a 
QSO absorption system at a level of (D/H)qqo « 1.9 — 2.5 x 10~ 4 (Songaila et 
al. 1994; Carswell et al. 1994). 

Comparison With Previous Bounds 

The bounds of primordial deuterium derived above are quantitatively simi- 
lar to previous bounds (YTSSO; WSSOK; ST92) but, there are some qualitative 
differences. One source of these differences is the new Geiss (1993) solar sys- 
tem abundances; another is the use of mass fractions (^2) rather than number 
ratios to hydrogen (y 2 ). For example, YTSSO had derived a version of equa- 
tion (8) which neglected X Q /Xp(^ 0.92) and related y 23 p (rather than y 2 p\ 
V2P < V23p) to 1723© and y 3Q : y 23P < 2/23© + (9 3 q _ !)l/3©- For the 2a upper 
bounds to 3 He and D + 3 He , YTSSO took: 7723© < 4.3 x 10" 5 , 773© < 1-9 x 10" 5 , 
which leads to: 10 5 t/ 2 3P < 2.4 + 1.9^0 . Then, for g 3Q > 1/4, YTSSO derive 
2/23P < 10 x 10 -5 which, corresponds to 7710 > 3 (YTSSO); actually, to two signifi- 
cant figures, 7710 > 2.8. WSSOK used the YTSSO constraint but, slightly different 
2a bounds to solar system D + 3 He and 3 He to derive: 10 5 y23P < 3.1 + 1.8^^,. For 
#3© > 1/4, WSSOK find y 23P < 10.3 x 10" 5 (WSSOK also derive the constraint 
y 23 pX P < 7.5 x 10~ 5 which, for Y P < 0.25, results in y 23P < 10 x 10~ 5 ); this, too, 
corresponds to 7710 > 2.8 (WSSOK). Finally, ST92 used the DSS stellar models for 
g 3 in the standard galactic evolution models of Tosi (1988) to identify the ranges of 
X 2 p and X 3 p allowed by the solar system and ISM data. For the 2a upper bounds 
to solar system abundances ST92 chose X 2Q < 4.5 x 10~ 5 and X 3Q < 3.6 x 10~ 5 ; 
for ISM deuterium they adopted X 2 jsm < 4.8 x 10~ 5 . In these models the effec- 
tive 3 He survival fraction is large (r 30 > 1/2, see the earlier footnote) and the 
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BBN constraint more restrictive: X 2 p < 9.0 x 10 5 ; 7710 > 3.7. Note that if the 
ST92 adopted abundances are used in eq. (10), 10 5 X° P < 4.2 + 2.4]?3q, which, for 
T3© > 1/4, implies X® P < 13.8 x 10~ 5 and r)% > 2.7; this is to be compared to 
X® P < 16 x 10~ 5 and 77^ > 2.5 using the Geiss (1993) abundances. Thus, although 
the Geiss abundances loosen the restrictions on rfl compared with YTSSO and 
WSSOK, it is seen that the neglect of primordial 3 He is significant and, the Geiss 
abundances in eq. (9) (with r 30 > 1/4), lead to the more restrictive constraint 
found above: X™ AX < 11 x 10" 5 ; ?7^ IN > 3.1. 



Summary and Implications 

The primordial abundances of D and 3 He cannot be derived from solar system 
or ISM data (with the exception that X 2 p > X 2Q , X 2 p > X21SM) without the 
intervention of a model to follow the galactic evolution of these nuclides. It has 
been the goal here to derive a "generic" analytic expression relating primordial and 
galactic abundances. Equation (9) incorporates the effects (and uncertainties!) of 
stellar and galactic evolution in one free parameter T 3 . For r 30 > 1/4 (DSS) and 
the Geiss (1993) solar system abundances it has been found, neglecting net stellar 
production of 3 He in low mass stars (Iben 1967, Rood 1972) that, X 2 p < X?f P AX < 
11 x 10" 5 . For standard BBN (WSSOK), this corresponds to y 2 p < 7.3 x 10" 5 and 
7710 > 3.1. This provides a slightly more restrictive lower bound to the universal 
abundance of nucleons (e.g. compare to WSSOK where 7710 > 2.8). For standard 
(N v = 3) BBN the lower bound to the predicted primordial mass fraction of 4 He 
increases from the WSSOK value to Ibbn > 0.241 ±0.001 as does the lower bound 
to primordial lithium: (Li/H) B BN > 1-2 x 10" 10 . For 7710 > 3.1, tt B > 0.045^ 2 so 
that, for H < 100 km s" 1 Mpc" 1 , tt B > 0.011. This strengthens, very slightly, 
the case for baryonic dark matter (Q B > Olum, see WSSOK). 

The upper bound to 77 from the lower bound to primordial D (X 2 p > 1.7 x 
10" 5 , 7710 < 9.0) is weaker than that from lithium [(Li/H) P < 2 x 10 10 , 7710 < 4.0] 
or from 4 He (Yp < 0.243, 7710 < 3.9). If this bound were saturated (7710 ~ 9.0), 
the primordial lithium abundance would be ~half solar and Yp m 0.253 ± 0.001 
(for N v = 3). 

In closing, it must be reemphasized that the analytic framework developed 
here is not a replacement for detailed evolution models. In the absence of such 
models, the choice of the 3 He survival fraction, r 3 , is little more than a guess. 
Rather, the potential value of this framework is that it provides a concise way to 
compare the results of different models; models which yield different constraints 
using the same observational data must differ in T 3 . Tracing the source(s) of these 
differences may help to choose among different evolution models and/or to decide 
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on the true uncertainties in the derived primordial abundances. 

We thank Daniele Galli for calling our attention to the issue of the correct 
meaning of the gs values published by DSS and adopted by ST92. This work was 
begun when G. S. was an Overseas Fellow at Churchill College, Cambridge and a 
Visiting Fellow at the Institute of Astronomy and he thanks them for hospitality. 
The work at OSU is supported by DOE grant DE-FG02-94ER-40823. 
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Figure Captions 



Figure la. The primordial deuterium mass fraction (X2P) as a function of the present 
nucleon-to-photon ratio rjiQ. The solid curve (.X^bbn) is the standard BBN 
prediction. The dotted curve (^max) an< ^ the dashed curve (^|max) are 
computed from equations (10) and (9) respectively using the 2a upper bounds 
to the abundances in eq. (13) for r 30 = 1/4. 

Figure lb. As in Figure la, for r 30 = 1/2. 

Figure 2. Upper and lower bounds to primordial deuterium. The solid curve (X 2 bbn) 
and the dashed curve (^max) are as m Figure 1. The dotted line (X^ lN ) 
is the 2<T lower bound to X 2 p from the ISM data (eq. 16). The restriction 
X 2Mm < -^2BBN < X%^ AX bounds 7710 : 3.1 < i] 10 < 9.0. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/astro-ph/9502067vl 



10 




^10 



This figure "figl-2.png" is available in "png" format from: 
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